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ABSTRACT 

The Doppler parameter distribution of Lya absorption is calculated for a set of differ¬ 
ent reionization histories. The differences in temperature between different reionization 
histories are as large as a factor three to four depending on the spectrum of the ion¬ 
izing sources and the redshift of helium reionization. These temperature differences 
result in observable differences in the Doppler parameter distribution. Best agreement 
with the observed Doppler parameter distribution between redshift two and four is 
found if hydrogen and helium are reionized simultaneously at or before redshift five 
with a quasar-like spectrum. 
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1 INTRODUCTION 

The physical conditions of the Intergalactic Medium (IGM) 
have been under discussion ever since its existence was pos¬ 
tulated. The absence of strong Lya absorption in the spectra 
of high-redshift objects showed that any homogeneous IGM 
must be either highly ionized or extremely tenuous (Gunn & 
Peterson 1965, Bahcall & Salpeter 1965). Nevertheless the 
IGM density and temperature have been poorly constrained 
for many years. Both a hot collisionally ionized IGM heated 
by energy input due to supernovae and a warm intergalac¬ 
tic medium heated by photoionization have been discussed 
extensively in the literature. Recent progress in numerical 
simulations including the relevant gas physics suggests that 
the low column density Lya forest is due mainly to density 
fluctuations of moderate amplitude in a warm photoionized 
phase of the IGM (Cen et al. 1994; Petitjean, Mricket & 
Kates 1995; Zhang, Anninos & Norman 1995; Hernquist et 
al. 1996; Miralda-Escude et al. 1996; but see also Bi, Borner 
& Chu 1992 for early analytical predictions). The baryon 
content in this warm phase must be similar to that predicted 
by cosmic nucleosynthesis to produce the observed mean ab¬ 
sorption in the spectra of high-redshift quasars. At redshifts 
larger than two most of all baryons are therefore probably 
in such a warm phase of the IGM and little room is left for 
an additional hot phase. (Rauch & Haehnelt 1995, Rauch et 
al. 1997, Weinberg et al. 1997). As pointed out by a number 
of authors the recombination time scale of the IGM becomes 
longer than the Hubble time at a redshift of about five. As 
a result, the temperature of the warm IGM has some mem¬ 
ory for when and how it was (re-)ionized (Miralda-Escude 


& Rees 1994, Meiksin 1994, Hui & Gnedin 1997). We have 
used a set of hydrodynamical SPH simulations for different 
plausible reionization histories to demonstrate that these do 
indeed result in significant and observable differences in the 
thermal history of the IGM. In this letter we investigate the 
observable consequences on the Doppler parameter distribu¬ 
tion of Lya absorption lines and speculate on the possibility 
to discriminate between different reionization histories. 


2 REIONIZATION AND THE THERMAL 
HISTORY OF A PHOTOIONIZED 
INTERGALACTIC MEDIUM 

We have calculated the thermal history of an IGM of pri¬ 
mordial composition within a cosmological context using 
SPH simulations performed with GRAPESPH (see Stein¬ 
metz 1996 for a detailed description of the numerical tech¬ 
niques and Rauch, Haehnelt & Steinmetz 1997 for a de¬ 
scription of their absorption properties). The code includes 
the relevant cooling and heating processes and follows self- 
consistently the non-equilibrium evolution of the baryonic 
species (H, H + , He, Hc + , He ++ , and e~). For the low-density 
regime relevant in this paper the thermal evolution of the 
IGM is rather simple (Miralda-Escude & Rees 1994, Hui & 
Gnedin 1997, Steinmetz & Haehnelt 1997). The dominant 
heat input is due to ionization of hydrogen and helium while 
the main cooling processes are Compton cooling at 2 <; 5 
and adiabatic cooling at lower redshift. At reionization the 
gas acquires an initial temperature 
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Figure 1. The left and middle panels show the absorption cross section-weighted specific intensity of the UV background at the Hi 
and Hen ionization edges ( J l u a^-l J cr-^) in units of erg cm -2 s —1 Hz -1 sr —1 for four different models. The right panel shows the 
corresponding thermal history of the IGM at the mean density. For a description of models A-D see section 2. 


Ti = 3^b / 4?r Iu [(^ “ hu Ri) CT Hi 

+ - ^Hei) a Rei + ( hv ~ ^Hen) CT Hen] 

/ J ^ I" [°Hi + CT Hei + °I1 cii] 7^7> 

(1) 

where ux and ax are ionization energies and ioniza¬ 
tion absorption cross sections and /„ is the specific intensity 
of the ionizing background. Afterwards the gas cools adia- 
batically within a Hubble time to the (redshift dependent) 
equilibrium temperature where photo-heating balances adia¬ 
batic cooling. The temperature of the gas therefore depends 
on the spectrum of the ionizing source and on the redshift 
when reionization occurs. 

We have investigated four different models to cover the 
uncertainty concerning the respective role of quasars and 
stars in the reionization of the Universe: 


Model A: 

Model B: 

Model C: 

Model D: 


Ionizing background as proposed by Haardt & 
Madau (1996). The spectrum is a power-law 
with a=1.5 processed by the intervening Lya 
absorption. 

Same as model A but the redshift evolution 
is stretched towards higher redshift (“Early 
QSO’s”). 

Same as model A but the redshift evolution 
is compressed towards lower redshift (“Late 
QSO’s”). A stellar component is added to the 
UV background which reionizes hydrogen at 
redshift z ^ 6 with a soft spectrum (power law 
with a=5). 

Same as C but the stellar component reionizes 
hydrogen at redshift a ~ 30. 


Figures la and lb show the evolution of the UV back¬ 


ground for the four scenarios at the Hi and Hen ionization 
edges respectively. At redshift three the Hi ionizing flux is 1 
to 3 x 1CV 22 erg cm -2 s^ 1 Hz^ 1 , slightly lower than estimates 
from the proximity effect but consistent with the mean flux 
decrement in QSO absorption spectra (Bechthold 1994; Gi- 
allongo et al. 1996; Cooke, Espey & Carswell 1997; Rauch et 
al. 1997). Model A was chosen by Haardt & Madau to repre¬ 
sent the UV background due to observed quasars. Model B 
mimics the existence of an as yet undetected population of 
quasar at redshifts beyond five. Models C and D address the 
possibility that the Haardt & Madau model overestimates 
the UV background due to quasars at redshifts larger than 
three. In the latter case the UV background would have to be 
dominated by a stellar contribution at these redshifts. Fig¬ 
ure lc shows the temperature evolution of the IGM for the 
mean density at a given redshift. There are striking differ¬ 
ences in the thermal history between the models. In models 
A and B hydrogen and helium are reionized almost simulta¬ 
neously with a rather hard AGN-like spectrum. This leads 
to an initial temperature of about 2 x 10 4 K. Thereafter the 
temperature approaches the redshift dependent equilibrium 
temperature. In scenario C and D the stellar component of 
the UV background is too soft to reionize helium at the 
same time as hydrogen. This results in an initial temper¬ 
ature of 7 x 10 3 K immediately after reionization - signif¬ 
icantly smaller than in model A and B. In Model C the 
late reionization of helium due to the QSO component of 
the UV background leads to a corresponding jump in the 
temperature while in model D Helium reionization occurs 
gradually before the QSO component of the UV background 
dominates. In the observationally relevant redshift range be¬ 
tween five and three the difference in temperature between 
the different models is as large as a factor 3 to 4. 
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Figure 2. The distribution of the Doppler parameter obtained by 
lines have column densities in the range 12.8 < IVjy ( < 16.0. 

3 SIMULATED AND OBSERVED DOPPLER 

PARAMETER 

In Figure 2 we show the distribution of the Doppler pa¬ 
rameter for the four scenarios at 0=2,3,4. The Doppler pa¬ 
rameters were obtained by fitting Voigt profiles to artificial 
spectra with the automatic line fitting program AUTOVP 
kindly provided by Rorneel Dave (Dave et al. 1997). At red- 
shift two the temperature in all four models is very similar 
and this is reflected by a nearly identical distribution of the 
Doppler parameters. At z > 3 the considerable differences in 
the temperature between the four models clearly carry over 
to the Doppler parameter distribution. In models C and D 
the distribution extends to smaller values of the Doppler 
parameter than in models A and B by about 5km s -1 . 

In the following we will present a preliminary compari¬ 
son between observed and simulated Doppler parameter. It 
should, however, be kept in mind that that the primary aim 
of our paper is to demonstrate that there are observable dif¬ 
ferences between different reionization histories. There are a 
number of uncertainties in any comparison with the observed 
data. The observed Doppler parameter distributions we dis¬ 
cuss have been obtained with different fit procedures. The 
spectrum of the UV background component due to star for¬ 
mation could be somewhat harder than the power law with 
a = 5 which we have assumed. Furthermore our simulations 
do not include any feedback effects due to star formation 
and may therefore underestimate the fraction of the IGM in 
a hot collisionally ionized phase. 

Figure 3 summarizes the comparison of observed and 
simulated Doppler parameter distribution at three different 
redshifts (Lu et al 1996, Kim et al 1997). The mean redshift 
and the column density range of the sample of lines identified 
with AUTOVP in the artificial spectra has been matched to 
those of the observed samples. At high redshift ((0)=3.3,3.7) 
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profile fitting artificial spectra for models A-D for z = 2,3, 4. The 


there is good agreement between the observed and simu¬ 
lated Doppler parameter distribution in the case of models 
A and B. The only discrepancy is a modest high-velocity tail 
which is present in the observed data but not in the simu¬ 
lated distribution. Models C and D clearly disagree with the 
observed distribution. They show a considerable fraction of 
lines which are cooler than the lower cut-off of 15 kms -1 . At 
redshift 2.3 there is a significant discrepancy between sim¬ 
ulated and observed distribution for all four models. The 
observed distribution shows a pronounced high-velocity tail 
which is not reproduced by any of the models. The existence 
of this high-velocity tail has been known for some time (e.g. 
Carswell et al. 1987) and appears to be real (Rauch et al. 
1992, Press & Rybicki 1993). 


4 DISCUSSION AND CONCLUSIONS 

We have demonstrated that the range of plausible reion¬ 
ization and thermal histories for the universe result in ob¬ 
servable differences in the Doppler parameter distribution of 
intermediate column density Lya absorption lines at 0 ~ 3 
to 5. 

The observed Doppler parameter distribution at 0 <; 3 
and the drop of the lower cut-off of the distribution from 20 
to 15 kms -1 between redshift 2.3 and 3.7 is well reproduced 
if the universe is reionized with an AGN-like spectrum where 
hydrogen and helium are reionized simultaneously. 

In the models where reionization occurs with a soft 
stellar spectrum roughly 10% of the lines at 0 ~ 3.5 have 
Doppler parameters between 10 and 15 kms -1 , well below 
the cut-off in the observed Doppler parameter distribution 
at this redshift (15 kms -1 ). If the universe were to be reion¬ 
ized by a stellar component of the UV background, either 
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Figure 3. Observed (dashed line) and simulated (solid line) Doppler parameter distributions for models A-D (from left to right). The 
mean redshift of the samples is (z)=2.3,3.3,3.7. The observed data are taken from Lu et al. (1997) and Kim et al. (1997). The histograms 
of the simulated Doppler parameter are renormalized in each panel such that the area under the curves is the same as that of the observed 
distribution. 


its spectrum would have to be hard enough to ionize helium 
and hydrogen at the same time, or the QSO contribution 
would have to become dominant before redshift four. 

None of our models can explain the observed high- 
velocity tail of the observed low-redshift sample {{z)= 2.3). 
About 15 % of all observed lines at this redshift have Doppler 
parameters of 40 to 60 kins -1 which would correspond to 
temperatures of 1 to 2 x 10 s K if the broadening were purely 
thermal. This is well above what can plausibly be reached 
by photoionization. This may indicate either, (i) that our 
simulations significantly underestimate the fraction of the 
IGM in a hot collisionally ionized phase because they do 
not include any feedback effects due to star formation or (ii) 
that there is a turbulent component to the broadening on 
scales below the resolution limit of the simulations which is 
more important at lower redshift. 

Our results suggest that a larger and homogenous sam¬ 


ple of lines, extending over a wider redshift range, will be an 
excellent tool to further constrain the reionization history. 
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